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In this work we propose a method for including strong isospin breaking in twisted mass lattice 
calculations, while preserving flavor identification. We utilize a partially quenched construction in 
which the sea quarks are given by the standard twisted mass lattice action while the valence quarks 
have an additional strong isospin breaking mass term. This construction allows for a practical 
use with existing twisted-mass gauge ensembles. Additionally, we construct the relevant partially 
quenched twisted mass chiral perturbation theory for both mesons and baryons to 0{m^,m q a,a 2 ). 
We provide explicit expressions for the pion, nucleon and delta masses, as well as the corresponding 
0^ ' mass splittings, and discuss the resulting errors from including the strong isospin breaking in the 

, valence sector only. Finally, we demonstrate how the application of this idea can be used, with mild 

■ approximations, to determine the values of both the up and down quark masses. 
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I. INTRODUCTION 



In the last few years, we have witnessed the development of twisted mass lattice QCD 0, [| grow from its initial 
quenched studies 0, 0, H, H, into a fully viable method for including two flavors of dynamical fermions d, [H, 
Hoi [TTI . [l2| . with more recent extensions to 2 + 1 + 1 flavors EH, E3- The twisted mass lattice action offers the 
promise of automatic O(a) improvement (where a denotes the lattice spacing), provided one can "tune to maximal 

£n- twist" m es d h os m. Further, the twisted mass term protects the Wilson-Dirac operator from acquiring 

zero or negative eigenvalues, as the symmetry properties of the twisted mass lattice action prevent additive mass 
renormalization, thus allowing for a cheaper numerical push to light quark masses. In a parallel effort, twisted mass 
chiral perturbation theory was developed for mesons [2ll . I22I l23l l24j as well as baryons [25[, both to aid in the 
understanding of the phase structure of the theory [26| and allow for a combined continuum-chiral extrapolation to 
' the physical point. For a nice review of twisted mass lattice QCD, see Refs. (27l. [28j. 

There are a few drawbacks however. First the action must contain an even number of fermion flavors. Second, the 
^■p ' action at finite lattice spacing breaks parity (although it has a combined parity and flavor transformation symmetry), 
! and thus the study of odd (even) parity states is complicated by their overlap with their lighter even (odd) parity 
• ■ counterparts. Third, the theory breaks the SU(2) flavor symmetry of li ght quarks down to U(l). As an extreme 
example, this allows for a mixing of the / = 2, 1% = and 1 = irir states |29j|. complicating or rendering impractical 

■ the calculation of the / = (I3 = 0) scattering channels. Lastly, the inclusion of strong isospin breaking effects 
I (m u 7^ md) leads to a non-perturbative mixing of the quark flavors 25|, thus complicating the computation of matrix 

• • , elements in which flavor identification is crucial. It is this last problem we address in this work. 
. £h ' The inclusion of strong isospin breaking in the twisted mass lattice action will allow for the computation of phe- 

, nomenologically interesting and important physical matrix elements. Most notably, it will allow for a lattice deter- 

■ initiation of the up and down quark masses [30(. Further, it will allow for a computation of the isospin violations 
5t i m the hadron spectrum, for example as was recently done for the neutron-proton mass splitting utilizing a partially 

quenched framework in which the sea quarks were degenerate [3lj . 

In this work, we propose a partially quenched extension to the standard twisted mass lattice action which introduces 
a strong isospin violating quark mass in the valence sector, with an average mass equal to the mass of the degenerate 
pair of twisted mass sea quarks. After detailing our proposed action in Sec. [TTJ we proceed to construct the relevant 
partially quenched twisted mass chiral perturbation theory for both mesons, to leading order (LO) and next-to-lcading 
order (NLO), as well as that for nucleons and deltas through next-to-next-to leading order (NNLO), Sec. IIIII In this 
section, we also discuss the vacuum alignment of the theory and the hairpin structure from the partial quenching. We 
then proceed to determine explicit expressions for the pion, Sec. IIII CI and nucleon and delta masses in Sec. IIII D[ as 
well as detail the mass splitting expressions. Importantly, we show in this partially quenched framework, the NLO 
corrections to the nucleon (and delta) masses exactly cancel. Therefore, the expansion of the nucleon mass splittings 
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should be as well behaved as for the pion mass itself. We then conclude in Sec. IIV1 showing how this method can be 
used to determine the values of the physical light quark masses. 



II. PROPOSED ACTION 



In order to prove the existence of a positive-definite quark determinant with a twisted-mass lattice action and 
mass non-degenerate flavors, it was shown the Pauli-matrix used for the isospin breaking must be different from that 
used with the twisted mass term (32|. However, as was explored in Ref. [2q . with a focus on baryons, because the 
discretization effects in hadronic quantities are expected to be approximately the same order of magnitude as the 
light quark effects, this leads to a non-pcrturbative mixing of the up and down flavors. For example, with a quark 
doublet q = (51,(72), an isospin breaking mass term implemented with T3 would lead to the identification q u = q\ 
and qd = qi while a twisted-mass term proportional to n would lead to q u = (qi + 92) and qd — (qi — q2)/V%- 
The non-perturbative mixing of these two bases can be resolved, for example, order by order in twisted-mass xPT by 
comparing to numerical results. 1 

The requirement of different Puali-matrices for the mass splitting and the twisted-mass term arises from the desire 
to have a positive definite fermion determinant. This restriction does not apply to the computation of fcrmion 
propagators, and naturally leads to the consideration of a partially quenched action. We begin with the target 
continuum theory. We propose an action for which the continuum limit is a partially quenched theory with non- 
degenerate valence quarks, with an average mass equal to degenerate pair of sea quarks. In terms of the valence-sea- 
ghost quark labels, the quark mass matrix is given by 



mQ 



(' 



\ 



\ 



5 J 



(1) 



and the continuum Euclidean Lagrangian is given by 

£f m Q =Q[p+ml-&rS]Q, 

with 

Q i^vah ^vali ^seai &seai ^ghostt dghost^) j 

7# = diag(l, -1,0, 0,1,-1). 



(2) 



(3) 
(4) 



We know the continuum Lagrangian of the standard implementation of the twisted mass lattice action can be written 



iu;T37 5 /2 



(5) 
(0) 



where q T = (u, d) and T3 is the standard Pauli-3 matrix. Writing this in more standard form 

C tm = q' [Jp+ m't + «/XT 3 7 5 ] q' . 
Therefore, the lattice action we propose can be determined from 

Z p t 2 = Qe iuT ^^ 2 [p+ ml - &rg] e iUT ^/ 2 Q , (7) 

with 



diag(l, -1,1, -1,1,-1). 



(8) 



This is not expected to be as much of an issue in the numerical implementation of the strange— charm doublet [T3I . H3 . This is because 
the mass splitting term which is used is much bigger than the typical discretization effects, which are typically on the order of the light 
quark masses. 
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At maximal twist, u> = tt/2, this Lagrangian is 

Cj£ = Q[lp+ ifxr^ 5 - i5Vvls] Q , (9) 



where 2 



7V = diag(l, 1,0, 0,1,1), (10) 



is the valence projector. This can be thought of a pair of non-degenerate Osterwalder-Seilcr valence quarks [331 ] . with 
masses equal to mos = \fii S\. 3 

From this point, there are two practical ways one can proceed, which are equivalent in the continuum limit. Option 
one is to add the following mass term during the calculation of the valence propagators; 

_ $ T v e iwr 3 t,s 75 = f -6e iujs for the up propagator , . 

3 1 +5e~ quJl5 for the down propagator ' ^ ' 

where u> is determined from the given twisted-mass lattice action with (5 = 0. The second option is to add the mass 
term, assuming perfect maximal twist has been achieved; 

.„ _ J —i5j5 for the up propagator , . 

' ~~ y —i5"/5 for the down propagator ' ^ ' 

While the second option may be simpler to implement numerically, the first option is slightly cleaner theoretically 
although we emphasize again both are equivalent in the continuum limit. We therefore proceed assuming option 
one (jlip has been implemented and will point out the differences if the second option were implemented. 

Before proceeding, we examine the symmetry properties of this new mass operator. The standard twisted mass 
action has an exact combined parity and flavor exchange symmetry, V F , 

Ua{x) — ► Ua(xp) , xp = (—x,t) 
U k {x) -» u£(x P ) , k=l, 2, 3 



(13) 

1p(X) -» lTi"/4l{j(Xp) 



The new mass term explicitly breaks this symmetry, 

- i5QVvl 5 Q — > +iSQVv1 5 Q ■ (14) 

Thus, for example, the twisted mass term will now receive additive corrections. However, this symmetry breaking 
is only present in the valence sector, and will not feed back into the gauge action. Furthermore, this breaking must 
vanish as <5 — > 0, and thus in non-parity violating matrix elements, the error must be of 0(S 2 ) or higher, including 
the additive correction to the twisted mass term, or the value of ^critical- 



III. PARTIALLY QUENCHED TWISTED MASS CHIRAL LAGRANGIAN 



To understand the ramifications of this lattice action, we study the low energy theory with partially quenched 
xPT [34l l35l l36l . l37l . [H, H^, 5l| . To begin, we must first construct the continuum Symanzik action at the quark 
level |42l. |43||. and then proceed to construct the resulting chiral Lagrangian |44j |. Our proposed lattice action is given 
by 4 

S = 5 ]T [ 7 *(V; + V„) - r V;V„] + m + i 15 ^t% s - S^e^'^ \iix) (15) 



2 Note that [r|,r| s ] = [t%,V v ] = [Tv,t£ s ] = 0. 

3 We thank Steve Sharpe for drawing our attention to this connection. 

4 The ghost quarks are not actually computed in the numerical simulations, but are a theoretical construct to account for the lack of 
closed valence quark loops. 
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Expanding about the continuum limit, the effective continuum Lagrangian for the quarks can be written (to LO in 
the lattice spacing) 



C = q 



>+ (m - StZ) e"" 1s +ac su , ia^F^ 



The renormalized quark mass is defined in the standard way with 

m = m cos(w) = Z m (m — m c ) /a , 
fi = m sin(o>) = Z^fi^/a . 



(16) 



(17) 
(18) 



Here, m c is the critical mass (defined up to an 0(a) shift). The twisted mass term, ifir^js, is protected from additive 
mass corrections from the symmetries of the action (up to the 0(S 2 ) corrections mentioned above). The twist angle 
is defined through the ratio tan(w) = fi/m'. Similarly, the valence isospin breaking mass term will also be protected 
from additive mass renormalization (to all orders), 



S = ZsSo/a. 



(19) 



There are additional operators at 0(a), but these are suppressed by additional powers of the quark masses and will 
not modify the construction of the chiral Lagrangian through NLO (24[. The resulting partially quenched meson 
Lagrangian is given through Q(rrvt, m q a, a 2 ) by 



C (PQ) = L str (fl EfyEt) - L- str (x' f E + £ V) - l[ PQ) [str (d^d^)] 2 - L (PQ) str (0 M E0 v £t) str (0 £$, E t) 

o o 

- L (PQ) str (a^Ea^Eta.Ea.Et) + L {PQ) str (a^Ea^Et) str ( X 'tS + EV) - L[ PQ) [str (x'ts + Et x ')] 2 
+ 4 PQ) str (a,Sa,Et ( x 't S + Et^)) - 4 PQ) [str (x*E - SV)] 2 - L (PQ) str ( X 'W^ + SV SV ) 



TV 



(PQ) 



str (a^Ea^Et) str (its + £U 



w, 



(PQ) 



str 



it£ + SU 



TV, 



(PQ), 



str (x /f E + Etx') str (its + £U) - TV 7 (PQ) str (x' f E - SV) str (it£ - EU 



TV. 



(•PQ) , 



str (x^SitE + E+x'sU 



TV, 



'(PQ) 



str f i f £ + £ti)l - W/' (PQ) str f i f Ei f E + E^SU 



(20) 



Here, wc have already absorbed the leading discretization effects into a redefinition of x [451 ] : 



X' = 2B ( m' + t>7f - 5T|e twr 3' 
A = 2W a = a 



2W a 



Further, we assume a power counting 



\J{rh' + a) 2 + (i 2 ~ <5~ a. 



(21) 



(22) 



The Lagrangian (|2T)|) is then the complete NLO meson Lagrangian relevant for our work. To extract relevant physics 
information, this Lagrangian must be matched to its non partially quenched counterpart. The unquenched Lagrangian 
is given by [24J 



/ 



£ = Y tr (^^Et) - ^-tr (x' f S + Et X ') - j [tr (d^d^)]* - |tr (a^Et) tr (d^tf) 
- ^ ^ (x' f E + Et x ')] 2 + | tr (a„E0„£t) tr ( x 't S + sV ) _ il [ tr _ gt^)] 2 



h 



+ TVtr (a^Ea^Et) tr ("itE + Eti) - W tr (x' f £ + £ f x') tr (i f E + E f i) - TV' 
- TV 7 tr (x' f £ - £ f x') tr (i f £ - 



tr ( i f E + E f i 



(23) 
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Here, the U are the SU(2) Gasser-Leutwyler coefficients [46|. The matching of the continuum NLO Lagrangian has 
been performed [47j | , and here we also match the partially quenched twisted mass Lagrangian; 

h=AL[ PQ) +2Li PQ) l 2 = ±L[ PQ) 

h + h = 164 PQ) + 8Li PQ) h = %L[ PQ) + AL[ PQ) 

l 7 = l64 PQ) +8Li PQ) 

w = wi PQ) + \w { b PQ) w = w ( / Q) + l -wi PQ) 

W? = w fQ) + X _ W iPQ) W > = w 'iPQ) + l _ w 'iPQ) (24) 



A. Vacuum Alignment 



To determine properties of the low energy theory, it is useful to first expand about the vacuum of the theory. In 
unquenched mass-degenerate twisted mass xPT, it has been shown that the vacuum is given at NLO by [24[ 



where lu = luq + e and 



exp(iw T3) 



S = (0|E|0) = exp(i(JT 3 ) 



(25) 



32 f Qi 

e(u) = — pro sin(w) ( W + 2W' cos(w) — 
/ \ m 

rh = \/ (m' + a) 2 + p? . 



(26) 



In the continuum limit, e — > as one would expect. 

In our proposed partially quenched theory, the vacuum has a slight additional perturbation. To find the vacuum, 
one can either minimize the potential energy (which is slightly more complicated for a partially quenched theories (48[) 
or one can require all single-pion vertices to vanish. Working through NLO, one can show that the Lagrangian, as 
written in Eq. (|20p . is rotated from the vacuum by the following angle 



with 



e"Vv) 



(27) 



where 



lu = cjo + e(w) 

«» 



m 2 — 5 2 
5m 



S 2 



and rh and e(w) are given by Eq. (|26p . and 

VT 



Jph = exp 



16asin(w)Fy 8 (PQ) 

16asin(w)W^ (PQ) 
P 



with 



(28) 



( Vu 


71+ 


<t>uj 


4>ui 


4>uu 


fad\ 




Vd 


4>dj 


4>dl 


4>du 


fad 


4>ju 


<j>jd 


>u 


<t>jl 


4> 3 u 


fad 


4>lu 


4>u 


<Pij 


m 


4>iu 


$l~d 


tyuu 


4>ud 


<j>uj 


<i>ui 




fad 


\<f>du 


fad 




<t>di 


fad 


fid J 



(29) 



Again, in the continuum limit, one finds the extra rotations of the vacuum vanish as expected. If we had used option 
2, Eq. (fT2"|) . then both e' and e" would have an additional O(a) shift from the tuning of lu. 
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B. Partially Quenched HairPin Interactions 

In the partially quenched Lagrangian, there are additional operators we have not written in Eq. (|20p . related to 
the singlet field <t> = str( ( f ) )/v r /Vp, which will ultimately be integrated out of the theory. These operators can 
be included to determine the structure of the hairpin interactions which give rise to the partially-quenched chiral 
logarithms [HI, H3, HI] . When calculating properties of pions, it is convenient to use the basis of fields (49l. [50j. 5 

o Vu - rj d _ i] u + ri d 

n= ^r> v = ^r- (30) 

One can show the propagators of these fields, which can be determined from the Appendix of Ref. (5lj | (this work 
considered a mixed action with twisted mass one of the possible types of sea quarks), 

Gn° = 2 , - , Qfj = , 2 , ZTo ■ (31) 

p -T m (p z + m) J 

With these propagators, one can determine the loop corrections to the pion masses and decay constants. The resulting 
loop corrections are given by 

where the loop function is defined as the regulated 4d Euclidean integral 

lI{ ™^L^tiT^> (33) 

which in the modified dimensional regularization of Ref. (4|| is 



This expression, with the known volume corrections to the one-loop pion mass formula [521 , [53j 

il[FV] - il[ooV] = V -r^K^nlm^L) , (35) 

4tH m^L \n\ 

can be used to trivially determine the volume corrections from the hairpin contribution (we have made use of the 
consistency at NLO of replacing the leading order pion mass with the full pion mass in these expressions, rh —> m 2 ,). 



Pion Masses 



Using this Lagrangian, we can now determine the pion masses. Working through 0(m 2 il m q a,a 2 ), and using the 



matching relations Eq. (I24|) . we find 



m ± = m 



mlo = mi ± 



(Air/)* 



+■ S^(A) 



A 2 



A 4 pq 
2(4tt/) 2 



32a cos(w) 

T 2 



m(2W -W) + 2W'acos{u] 



4<f 



&4W'a 2 sin z {uj) 
P 



(36) 



Keeping the mass of the singlet explicit (tno), these two fields are given at leading order in <5 by 7r° = — ~ d + 2 T? "^ ld and 
fj _ Vu+Vd _ 5 i" as t ne isospin breaking mass term allows a coupling between the ir° and rf fields. These extra terms need not 
be included in the propagators as they are included as interactions in the partially quenched chiral Lagrangian. 
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The first non-standard term in this expression arises from the hairpin interactions, and is the remnant of the partially 
quenched chiral Log. Because the fj propagator has extra suppression, Eq. ([3~Tj) relative to the standard version of this 
partially quenched propagator [49[, (this happens because rh va i = m sea ), the enhanced chiral logarithm has become 
simply a constant. To clarify which contributions arc physical, and which are partially quenched artifacts, we have 
introduced the term 



A 2 

^PQ 



S, 



(37) 



which we shall use in the remaining mass expressions. For 6 — ► these expressions reduce to those of standard 
twisted mass xPT. Therefore, with multiple values of the isospin breaking mass term <5, one can determine l-j from 
the charged-neutral pion mass splitting. At this order, the pion decay constants do not receive any corrections and 
are given by the standard form 



U = f 



1 - 



2m / rh 

(MP "I A 2 " 



2m I _, 1 , l6Wa 
-j^-kiA) +COSW— — 



(38) 



D. Baryons 

One can also include baryons in twisted mass xPT [25[ , using an extension of the heavy baryon chiral Lagrangian 
formulated by Jenkins and Manohar [54l , |55| . For our work, we will need the two-flavor partially quenched baryon 
Lagrangian, which was first developed in Ref. [4l[ and later extended to NNLO in Ref. (56|. Here we use the 
normalization conventions of Ref. [57] , for which the twisted mass baryon chiral Lagrangian is given by 



R {PQ) JPQ) APQ) 

= (Bv • DB) + (BBM + ) + g^y (BM+B) + ^ (BB) tr(M + ) + ^ (BB) tr(W + ) 

(PQ) —(PQ) —{PQ) 

(7> • DT^) - A (7 M T M ) + jf— (T,M + T,) — (T M T M ) tr(M+) - §^ (T„T M ) tr(W + ) 

(4tt/) (4tt/) (4tt/) 

+ 2a (PQ) (BS^BA^) + 2(3 {PQ) (BS^A^B) + 2H {PQ) (T'S^A^) + j|c [(T U A U B) + {BA V T V )\. 

(39) 

In this Lagrangian, ( ) denote the graded summation of flavor indices, first defined in Ref. :39j. The spurions arc 
defined as 



(40) 



with x' and A defined in Eq. ([21]) . Here, £ 2 — £ is needed for the inclusion of the baryon fields in the chiral Lagrangian. 
The axial field is defined as 



A, = -(^~eo^) , 



(41) 



and is a spin operator [54.l55|. As with the mesons, we must match this Lagrangian to the unquenched one, given 
by 



C = NV ■ DN + (ST) NM+N + (ST) NNt < M ^ + J^f) NNtr(W + ) 

+ (7> • D Tfj_) + A (T^) + (T^M+T^) + (T^) tr(M+) + (T^) tr(W+) 

+ 2 g A NS ■ A N - 2 g A A T^S ■ AT„ + g AN A f e n' N k + h - c ] • ( 42 ) 
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Performing the matching, one finds 



_ 2 (PQ) 1 R (PQ) 
a M - 3 a A/ - 3 Pm i 

ctm — er M + -a M + -p M , 

7a/ = 1 M > ct a/ , 



<M = §a<"» - \(i {PQ) , .91 = \cf™ + iffQ), 

.9 A A = W , ,9aw = -C, (43) 

and 

°w = <J W , °w = <J W , (44) 

for the discretization operators. Using the partially quenched Lagrangian, performing the matching with Eqs. (|43|) 
and (|44p . and working consistently to NLO we find the nucleon masses are given by 

„2 o„2 



m l ( a M . r ( K \\ 3 8 9AN t( a K \ 



Ol M 

(47T/ ff ) 2 + (4^)V 2 + <7 ^ (4^)2 m - 3(4 7 r/ 7r ) : 
acos(u)(j w 3nA 4 PQ (g A + gi) 2 

+ (471^) + 877^(47^)2 

J_°M , "4 /«M , r , A A 37r gj 3 Bg 2 AN 



acos(a>)o-iy 37rA^ Q (g A +.gi) 2 

+ (4ttA) + 877^(47^)2 
with 

a an /a2 2 ■ x3/2, / A+ N /A 2 -m 2 + ie\ 3 A 2l /m 2 \ a3i /4A 
^(m, A, A) = (A - m 2 + -) 3 / 2 In - _ ^= _ = = = - , Am ln ( A^ j " * 



(45) 



Here, we see that the NLO contributions exactly cancel in m n — m p , thus rendering the expansion of the mass splitting 
on the same footing as the pion mass expansion. One can see the last term in these mass expressions is proportional 
to the coupling of the nucleons to the singlet field, being proportional to (gA +<7i). These terms are remnants of our 
partially quenched theory and would vanish if the sea quarks had an isospin breaking mass term equal to that of the 
valence quarks. 

Similarly, one can determine the delta mass expressions. One should caution that due to the strong coupling to 
the tv — N system, the deltas, at lighter pion masses, have significantly larger volume dependence than the nucleons 
or pions [H, HH, [6(| ■ Neglecting these issues for this work, the delta mass extrapolation formulae are given by 



m- , A , S C A7A/ . mi pM . - r , . -A 25g^ A 3 2g 2 AN 
m A = M + A + — — + I — + cr M A) J - — ^ m% - a " %, -A, /i) 

(47r/ w ) 6 (47r/ 7r ) V 2 J 27(47r/ 7r )2 3(4tt/) 2 



acos(w)cr W ' 57r ^PQ5AA 



where the coefficients ca arc given in Table |T] and 



(46) 



^Ym -A A) - I A ' A) ~ 2 ^ (A2 ~ m2)3/2 ' m < 1 A| (47) 

S(m, A ' A J - | -F{m, A, A) - 27r(m 2 - A 2 ) 3 / 2 ; m > |A| ' 10 

In the limit 8 — > 0, these expressions reduce to those in Ref. (25|. Similar to the nucleons, the NLO contributions 
exactly cancel from the mass splittings, and the last term in this expression arises from the partially quenched 
construction. Also, at this order, one sees the delta masses obey an equal spacing rule, which is violated at NNLO by 
one operator [56| . It is precisely the imaginary piece of this ^"-function, which in finite Euclidean volume gives rise 
to the power-law dependence of the delta masses [59l . [60j . 
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TABLE I: Coefficients of the delta mass corrections arising from the LO term proportional to the strong isospin breaking and 
the NNLO discretization effect. 





A++ 


A+ 


A 


A" 


CA 


-3 


-1 


1 


3 




1 


-1/3 


-1/3 


1 



1. NNLO 



For the baryons, we can carry this prescription to NNLO (which is 0(m 2 ,m q a,a )). The complete set of partially 
quenched operators relevant for the nuclcon and delta masses at 0{m 2 ) was determined in Ref. [56[. The complete 
set of new twisted mass heavy baryon x?T operators at 0(m q a,a 2 ) was determined in Ref. [2^|. The extension 
of the twisted mass operators to the partially quenched theory is straightforward, however due to the cumbersome 
length, we do not detail them here. Rather we present the results of the mass corrections after the matching to the 
unquenched theory has been performed, and we only provide the expressions at maximal twist. For the nuclcons, 
there are 9 relevant operators in the partially quenched theory, which reduce to four in the unquenched theory, while 
for the deltas there are 7 relevant operators in both, 

C M = J^Jf^i 1 ^ M + N + bl I NNtx{M 2 + )+bl I N M+N tr{M+) + b* 1 NN [tr(X+)] 2 

+ if ff(M + M + y'T^ jk + if ff(M+y' (A4+)/ + *f T M T M tr(Mi) 



f (f^M+Tj ti(M+) + if f^[tr{M + )] 2 



+ b?~ NN tr(W_ W-) + ti (f^triW-W-) +t^-f^ i (W-) i i T^ yk . (48) 



As discussed in Refs. [1, the symmetries of the twisted mass lattice action prevent the twisted mass term from 
splitting the nucleon masses, which is reflected in the low energy chiral Lagrangian, while the delta-multiplet splits 
into two degenerate pairs [25| . The corrections to the delta masses from the twisted mass discretization, at maximal 
twist, are given by 

~ 2 t W ~ 

where c~a is given in Table [IJ 

With our particular partially quenched construction, the nucleon and delta masses will have an error at this order 
from one operator each, proportional to the terms with tr(.M?_). In the full theory with isospin breaking in both the 
sea and valence sector, the masses would receive mass corrections 



5m 7v 



2(47T/ ff ) 



3 



2(4tt/ w )3 

while in our partially quenched theory, these mass corrections become 



5m 



N 



2(47T/ ff )3 ' 



However, in the mass splittings, these contributions vanish leaving the baryon mass splittings free of strong isospin 
breaking errors at this order. At maximal twist, we then find the nucleon mass splittings is given through NNLO by 



10 



TABLE II: Coefficients of the delta equal mass splitting violation. 



ct. 2 



m A+ — m A++ 
-1 



m A+ 



thc expression 



(47T/.) 



OLM 



(4tt/ w ) : 



A I 



76 ' (4tt/ w )2 



45AjV ( g7M ™ OL M 



(47T/0 2 

a MAf> Q ^ o 3 a 



^"YmsUn - 4a M ( 5 ^ + .gljv) - a M (6 5 i + 1) In ( ^£ 



(52) 



where 



J(m,A,A) = 2Av / A2 



zeln 



VA 2 



Va 2 



to + le 



+ m In 



m 
A 2 



2A 2 In 



\ mr 



(53) 



Note, the nucleon mass splitting is free of discretization errors at this order. The last term in this mass splitting 
expression is not physical, having arisen from our partially quenched construction with Apg = 5, Eq. (|37|) . 

Of phenomenological interest, the delta mass operator with coefficient t^ 1 leads to a breaking of the delta equal 
mass spacing (56|. Taking the delta mass splittings between neighboring members of the multiplet, one finds (at 
maximal twist) 



iiia - mA' = 



f lM 



tf + tf J(m^, -A, A) g\ N {2 lM - a M ) 



(4tt/ w ) t 3 (AnUY 



(4tt/») 



M 
2 "-2 



52 7m 2 
81 3 5aa 



7mA 4 



18m2(47r/ T ) 2 



12 - 5<? AA 



7a/ 
3 

,w_ 

c t2 i 2 



(47T/ 7r ) ; 
?AA 



10^ 



27 



In 



A 2 



(47T/.) 



(54) 



where 



J(m,-A,A) 



J(m, A, A) + 4ittAVA 2 - m 2 , to < |A| 
J(m, A, A) - 4ttAVto 2 - A 2 , to > |A| 



(55) 



and Ct 2 is given in Table ITT1 In the limit 5 — > 0, the delta mass splitting reduces to that in Ref. [If 



TO A — TUA' 



ct 2 t 2 



(4tt/.)S 



(56) 



IV. DISCUSSION 



In this work, we have proposed a partially quenched extension of the twisted mass lattice action which includes a 
strong isospin breaking mass term in the valence sector, and preserves flavor identification. Further, we developed 
the corresponding partially quenched twisted mass chiral perturbation theory relevant for our proposed lattice action. 
While this proposal does not completely incorporate the effects of the strong isospin breaking, we have demonstrated, 
using the partially quenched theory, that the error in this approximation is sub-leading. This method extends to 
the hypcron spectrum in a straightforward manner. The relevant partially quenched twisted mass xPT could be 
created by adding the isospin breaking mass term to the valence quarks in Refs. [40l l6ll . [63 ] , and then matched to 
the unquenched theory [63| . Alternatively, one could develop the partially quenched version of the two-flavor xPT for 



11 



hypcrons [57j . One application of this work would be the calculation of a isospin violations in the hadron spectrum 
utilizing twisted mass LQCD, for which these errors are even further suppressed. Combined with the physical point, 
this would allow for a lattice determination of the physical light quark mass splitting. 26 = md — m u , a stated challenge 
in Ref. [3(j. 

To contribute to this proposed endeavor, we have determined the relevant chiral extrapolation formula for the 
mass splittings of the pions, nucleons and deltas (as well as expressions for the masses themselves). The leading 
non- vanishing correction to the pion mass splitting is 



a 2 4<5 2 MW >g2 sin 2 M 
m n o - m v± = —l 7 -2 . (57) 

J 7T J 



At a fixed lattice spacing, with multiple values of the isospin breaking mass parameter, one can determine the Gasser- 
Leutwyler coefficient, lj. Using the physical values of the charged and neutral pion masses, after correcting for 
electromagnetic interactions, this would allow for a lattice determination of the up and down quark masses. We 
have also provided expressions for the nucleon and delta mass splittings, and shown that through 0(to 2 , <5 2 ), these 
mass splittings are free from errors arising from the degeneracy of the sea quarks. Scaling the masses by one can 
perform an extrapolation of these mass splittings which requires no scale setting. The leading-order splittings (the 
NLO corrections exactly cancel and the full NNLO expressions are in Eqs. (|52|) and (|54[) respectively), 



m n — m p a^j 5 



O(Sml) 



0(Sml). (58) 



After correcting for electromagnetic effects, these baryon mass splittings, combined with the physical mass splittings, 
could also be used to determine the physical mass splitting of the light quarks (perhaps we should caution that for 
present lattice simulations, the nucleon mass behaves linearly in [64l. l65j| . a trend which is not presently understood 
theoretically, and may complicate the baryon extrapolation). Either way, one must determine the value of 5, both 
for a given set of lattice input parameters and at the physical point. A direct lattice calculation of this quantity 
requires the determination of the renormalization coefficient, Zg. However, one can make a good approximation of 
this parameter by taking advantage of the fact both the twisted mass parameter as well as the isospin breaking 
mass parameter are protected from additive mass renormalization (actually, the explicit breaking of the parity-flavor 
symmetry by the proposed isospin breaking mass parameter gives rise to an 0{8 2 ) correction to /it, but this is expected 
to be perturbative). From Eq. (f2"Tj) . one has 

I = (59) 
M z n Mo 

and at or near maximal twist, to a very good approximation, m 2 ~ jl. If one further approximate Zg ~ Z^, then a 
conservative estimate can be made for the isospin breaking mass parameter 

S^^ml, (60) 
Mo 

where So and (Iq are the bare mass parameters of the action, Eq. (|15p . A determination of the physical value of this 
parameter would allow a precise determination (supplemented by the approximations we have made) of the up-down 
quark mass splitting, 



5 m d - m t 



ft m d + m u { \( 47r A) 



l + O 77TW • (61) 
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